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A two dimensional array of Josephson junctions in a magnetic field is considered. It is 
shown that the dynamics of the vortices in the array resembles that of electrons on a two- 
dimensional lattice put in a magnetic field perpendicular to the lattice. Under appropriate 
conditions, this resemblance results in the formation of a quantum Hall fluid of vortices. 



The bosonic nature of vortices and their long range logarithmic interaction make some of 
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the properties of the vortices' quantum Hall fluid different from those of the electronic one. 
Some of these differences are studied in detail. Finally, it is shown that a quantum Hall 
fluid of vortices manifests itself in a quantized Hall electronic transport in the array. 
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1. Introduction 

This paper discusses a quantized Hall effect (QHE) state of vortices in a two dimensional 
(2D) array of Josephson junctions. Motivated by the analogy between Magnus force acting 
on a vortex moving in a two-dimensional ideal fluid and Lorenz force acting on a charge 
in a magnetic field, we study the transport of vortices in a Josephson junction array. In 
particular, we focus on the case in which the charging energy of the array is minimized when 
the number of Cooper pairs on each element of the array is not an integer. We find that 
for a certain range of parameters the vortices are expected to form a quantum Hall fluid, 
and the resistivity of the array is expected to show a QHE behavior. While some of the 
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properties of the quantum Hall fluid formed by the vortices are similar to those of the well 
known Laughlin fluid, formed by electrons in QHE systems, we find that the logarithmic 
interaction between the vortices leads to interesting modifications of other properties. 

The paper is organized in the following way: in Section (2) we review the classical and 
quantum mechanical analogies between the 2D dynamics of charged particles under the 
effect of a magnetic field and that of vortices in a 2D fluid. These analogies, arising from 
the analogy between Magnus and Lorenz forces, motivate the introduction of the system we 
analyze - a Josephson junction array in a magnetic field, and the study of a quantized Hall 
effect in that system. Section (2) is concluded with a precise formulation of the problem 
to be studied. In Section (3) we analyze the transport of vortices in this array. We show 
that the dynamics of the vortices can be mapped on that of charged particles under the 
effect of a magnetic field, a lattice-induced periodic potential and a mutual interaction. The 
mutual interaction is composed of a logarithmic " static" part as well as a velocity-dependent 
short ranged part. In Section (4) we analyze the quantized Hall effect associated with the 
transport of the vortices, and its observable consequences. In particular, we study the unique 
features of the QHE for logarithmically interacting particles. Conclusions are presented in 
Section (5). 

The possibility of Quantum Hall phenomena in Josephson junction arrays was recently 
discussed in two other works, one of Odintsov and Nazarov and the other of Choi @. 
The regime of parameters we consider is different from the ones considered by these authors. 
We comment briefly on this difference and its implications in section (2). 

2. Transport of vortices in a Josephson junction array — introduc- 
tion and motivation 

The classical dynamics of vortices in two-dimensional ideal fluids is well known to re- 
semble that of charged particles under the effect of a strong magnetic field || . An electron 
in a magnetic field is subject to Lorenz force, while a vortex in an ideal fluid is subject to 
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Magnus force. Both forces are proportional and perpendicular to the velocity. Two electrons 
in a strong magnetic field encircle each other, and so do two vortices in an ideal fluid. The 
dynamics of both are well approximated by an Hamiltonian that includes only a potential 
energy V(x, y), where x, y are the planar coordinates, and for which x and y are canonically 
conjugate. This approximation is known as the guiding center approximation for the elec- 
tronic problem, and as Eulerian dynamics for the vortex problem. This close resemblance 
naturally raises the possibility of analogies between transport phenomena of electrons in a 
magnetic field and those of vortices in ideal fluids. 

A vortex in a fluid can be viewed as an excitation in which each fluid particle is given 
an angular momentum / relative to the vortex center |4]]. Consequently, the velocity field 
v(r) of the fluid satisfies J r v-dl = — , where m is the mass of a fluid particle, and T is 
a curve that encloses the vortex center. When the vortex center moves with a velocity u, 
and the fluid is at rest far away from the vortex center, the vortex center is subject to a 
Magnus force, given by F Magnus = lulnu x z, where n is the number density of the fluid far 
away from the vortex core. Being both proportional and perpendicular to the velocity of 
the vortex center, Magnus force obviously resembles the Lorenz force acting on an electron 
moving in the x — y plane under the effect of a magnetic field Bz. This Lorenz force is 
given by FL orenz = —u x z, where u is the velocity of the electron. Thus, the role played 
by the product - c B in the latter case is played by the product 2nln in the former. While 
the fluid density plays a role analogous to that of a magnetic field, a fluid current plays a 
role analogous to that of an electric field. To see that, note that in a frame of reference in 
which the electron is at rest, the force it is subject to looks as if it arises from an electric 
field, given by — u x z. Similarly, in a frame of reference in which the vortex center is at 
rest, the force it is subject to seems to arise from the motion of the fluid. Since the fluid 
current density is J = nu, in the vortex rest frame the Magnus force is F Magnus = 2irlJ x z, 
and J x z plays a role analogous to that of an electric field. Thus, while the fluid density 
affects the vortex dynamics in the same way a magnetic field affects electronic dynamics, 



3 



the fluid current plays the role of an electric field. Maxwell's equation Vxi? + 4j| = 0is 
then analogous to the continuity equation in the fluid V • J + ^ = 0. || f| 

Quantum mechanics introduces two new ingredients to the analogies discussed above. 
The first is the quantum of angular momentum /, given by U (or alternatively, the quantum 
of vorticity, — ) @. The second is the quantization of the magnetic flux, the integral of the 
magnetic field over area. This quantization is most clearly seen through the Aharonov-Bohm 
effect ||: the Aharonov-Bohm phase shift accumulated by an electron traversing a closed 
path in a magnetic field is 2n times the number of flux quanta it encircles. Combining these 
two ingredients together, one should expect a quantization associated with the integral of 
the number density over area, i.e., with the number of particles. This quantization should 
manifest itself in the phase accumulated by a vortex carrying a single quantum of angular 
momentum, h, when it traverses a closed path in a fluid. Indeed, as shown first by Arovas, 
Schrieffer and Wilczek |§, such a vortex does accumulate a geometric (Berry) phase jlOj , 
and this phase is 2n times the number of fluid particles it encircles. The analog of a flux 
quotum is then a single fluid particle [llj . Note that the analogy between vortex dynamics 
in a fluid and electron dynamics in a magnetic field does not depend on the fluid being 
charged, and is valid for neutral fluids as well. 

Quantum transport of 2D electrons in a magnetic field crucially depends on the electronic 
filling factor, the ratio between the density of conduction electrons and the density of flux 
quanta. For a very low filling factor, (<^C 1), electrons are expected to form a Wigner lattice. 
At higher filling factors, the quantized Hall effect takes place ||12|| . Similarly, we expect 
transport of vortices in a fluid to depend on a vortex "filling factor", the ratio between the 
density of vortices and the density of fluid particles. However, in continious two dimensional 
fluids this ratio is usually much smaller than one, and the vortices indeed form an Abrikosov 
lattice. 

How can the vortices "filling factor" be made larger? In this work we make the vortex 
filling factor larger by considering a lattice structure. As is well known, all properties of 
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electrons on a lattice are invariant to the addition of a magnetic flux quantum to a lattice 
plaquette. Similarly, when we consider vortex transport on a lattice-structured fluid, we 
find all properties to be invariant to the addition of a single fluid particle to a lattice site. 
It is this periodicity that allows us to make the effective filling factor much larger than the 
ratio between the density of vortices and the density of fluid particles. 

Based on the foregoing general considerations, we study in this paper a Josephson junc- 
tion array in a magnetic field. Josephson junction arrays were extensively studied in recent 
years |L3| |TJ| . The array we consider is composed of identical small super-conducting dots 



coupled by a nearest-neighbors Josephson coupling Ej, and by a capacitance matrix C. 
For definitness, we consider a square lattice of the superconducting dots. Generalizations 
to other lattices are straight forward. A perpendicular magnetic field induces vortices in 
the configuration of the superconducting phase. We denote the average number of vortices 
per lattice plaquette by n v . Each of the dots carries a dynamical number of Cooper pairs, 
denoted by rij (for the i'th dot), as well as positive background charges. The charging energy 
of the array is minimized for a certain set of values of rij, which we denote by n Xti . (In our 
notation, charge is always expressed in units of 2e, dimensionless.) While the 

rij's are operators with integer eigenvalues, real parameters, that are closely related 

to the chemical potential of the dots. In this work we consider the case in which for all sites 



n x,i — n x- The Hamiltonian describing the array is [14] , 

H = ^£fa - n x )Cr\n 3 - n x ) + £j]T(l - cos(& - ^ - j° A ■ dl)) (1) 

1 ij (ij) Ji 

where denotes a sum over nearest neighbors, is the number of Cooper pairs on the 
i'th dot, (pi is the phase of the superconducting order parameter on the z'th dot, A is the 
externally put vector potential and the integral is taken between the sites i and j. A factor 
of — is understood to be absorbed in A. The matrix C~ x is the inverse of the capacitance 
matrix C. Generally, the matrix C includes elements coupling a dot to its nearest neighbors, 
to the substrate and to neighbors further away. The matrix elements of both C and C^ 1 
are a function of the distance between the sites i and j. For short distances the electrostatic 



energy is determined by nearest neighbors capacitance only, and all other capacitances can 
be ignored. The ij matrix element of C _1 is then ■^ ]L - log |r$ — r 7 -i, where C nn is the nearest 



neighbors capacitance For large distances, the electrostatic interaction depends also 

on capacitance to the substrate and capacitance to neighbors further away. The inverse 
capacitance matrix then decays with the distance. Throughout most of our discussion we 
assume that the size of the array is small enough such that the charging energy is determined 
by nearest neighbors capacitance only. Then, the charging energy involves one energy scale, 
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Eq = 2§ — . The effect of other capacitances is briefly discussed in section (3). 

Since our main interest in this study is focused on transport phenomena of vortices, we 
constrain ourselves to arrays in which Ej>Ec- In that regime of parameters vortices are 
mobile enough not to be trapped within plaquettes, but their rest energy is large enough such 
that quantum fluctuations of vortex-antivortex pair production can be neglected. Arrays 



in which Ej>Ec were studied experimentally by van der Zant et.al. ||15|| , and were found 
to show a magnetic field tuned transition from an almost super conducting state to an 
almost insulating state. At weak magnetic fields the density of vortices is low, and their 
ground state is the Abrikosov lattice. The array is then super-conducting. The transition 
to the insulating state, at a critical value of the magnetic field, is interpreted as caused by 
a transition of the vortices from a lattice phase to a correlated super-fluid-like phase [Tl| 



As mentioned in section (1), QHE phenomena in Josephson junction arrays were dis- 
cussed in two recent preprints. The first, by Odintsov and Nazarov focuses on the 
regime E c 3> Ej, and discusses a quantum Hall fluid of Cooper-pairs. The second, by Choi 
0, focuses on the regime Ej 3> Ec, and discusses a quantum Hall fluid of vortices. The 
quantum fluid we discuss in this paper has some similarity to the one discussed by Choi. 
However, the difference in the regime discussed, as well as our detailed study of the effect 
of the logarithmic vortex-vortex interaction, make some of our conclusions different from 
those of Choi. 

Due to the lattice structure of the array, the spectrum and eigenstates of the Hamiltonian 



(H) are manifestly periodic with respect to n x , with the period of one Cooper-pair (n x = 1). 
This periodicity is similar to the periodicity of the spectrum of electrons on a lattice with 
respect to the addition of one flux quantum per plaquette. Thus, although the ratio between 
the density of vortices and the density of charges in the system is very small, the physically 
meaningfull ratio is the ratio of n v to (n x — [n x ]) (where [n x ] is the largest integer smaller than 
n x ), and this ratio is not necessarily small. Following this observation, we limit ourselves 
from now on to the case < n x < 1. 

Having described in detail the Josephson junction array to be considered, we conclude 
this section by formulating precisely the question to be studied, namely, how do the physical 
properties of the array depend on the ratio between the vortex density n v and the charge 
density n x ? We start our examination of that question by deriving an effective action for 
the vortices in the array. 

3. The effective action for the vortices 

The Hamiltonian ([I]) describes the Josephson junctions array in terms of the sets of 
variables {rz,}, {(pi}. In this section we derive an equivalent description of the array in terms 
of the vortex density p vor , the vortex current J vor and gauge fields the vortices interact 
with. Our goals in attempting to derive this description are three-fold. The first goal is 
to verify the validity of our assertion that vortices are subject to Magnus force, and that 
n x plays a role analogous to that of a magnetic field in electronic dynamics. The second 
goal is to study the mutual interactions between vortices. The third goal is to estimate the 
mass of the vortices. The first two goals are relatively easy to achieve. Estimating the mass 
of the vortex, however, turns out to be a harder task, which we are able to handle only 
approximately. 

The effective action for vortices in a Josephson junction array was first discussed by 
Eckern and Schmid, who considered the Hamiltonian (|l|), with n x = 0. More generally, the 
effective action for singularities in the phase configuration of a complex field was discussed 
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in various other contexts in physics. A particularly convenient method for the derivation 
of such an action is the "duality transformation", developed and used by Jose et.al. ||16|| , 



Berezhinskii Peskin jnj, Fisher and Lee |19| and others. This method was applied to 
analyze the motion of vortices in Josephson junction arrays (again, for the case n x = 0) by 
Fazio, Geigenmuller and Schon PD|| . 

In our derivation of the effective action, we follow Fazio, Geigenmuller and Schon |2(| by 
applying the duality transformation to obtain an effective action for vortices on a lattice. 
The action resulting from the duality transformation (Eq. (§) below) describes the vortices 
as bosons on a lattice interacting with an externally put vector potential, as well as with a 
dynamical vector potential. The externally put vector potential, which we denote by K ext , 
satisfies V x JC ext = 27rhn x . The interaction with the dynamical vector potential mediates a 
vortex-vortex interaction, composed of two parts. The first part is the familiar logarithmic 
interaction. Its strength is proportional to the Josephson energy Ej. The second part, 
induced mostly by the charging energy of the array, is a short ranged velocity-velocity 
interaction. The latter makes the vortices massive, since it includes a self interaction term, 
quadratic in the vortex velocity. However, the mass defined by this interaction is a "bare 
mass" , that does not take into account the periodic potential exerted on the vortices by the 
lattice. Generally speaking, the periodic potential changes the bare mass into an effective 
band mass. In an attempt to estimate the band mass we write the continuum limit of 
the vortices action. In the continuum language, vortices are massive particles interacting 
with an external vector potential, with a periodic lattice potential and with one another. 
Our analysis of this rather complicated dynamics follows the way the dynamics of electrons 
on a lattice is analyzed. We start by neglecting vortex-vortex interactions. We are then 
faced with a single particle problem, in which a massive vortex interacts with a static vector 
potential JC ext , and with a periodic lattice potential. This problem is identical to the problem 
of an electron under the effect of a uniform magnetic field and a lattice periodic potential, 
whose solution is well known. When n x C 1 the effect of the periodic potential can be 
accounted for by changing the "bare mass" to an effective mass. We limit ourselves to this 
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case, and estimate the resulting effective mass. Then, we incorporate the vortex-vortex 
interactions back into the action. 

Before turning into the details of the derivation sketched in the last paragraph, we pause 
to define a notation. We denote 3-vectors by bold-faced letters, and their two spatial 
components by vector arrows. The electromagnetic potential is then A = (A , A x , A y ) = 
(Aq, A). We number array sites by a subscript i. The bond connecting a site i to its neighbor 
on the right side is denoted by the subscript i,x. Similarly, the bond connecting the i'th 
site to the site above it is denoted by the subscript i,y. The difference operator A, a 
discretized version of V, is defined accordingly. When operating on a scalar 0, for example, 
the x-component of Aj is (pj — (pi where j is the neighbor to the right side of i. 

Our derivation of the effective action starts by considering the partition function 

Z = tie^ H = J D{n t } J D{^} e -is(K(')}.W.W}) (2) 



where the action {&(*)}) = fg dt i hn^t^t) - H({m{t)}, {&(*)}) 



and the 



Hamiltonian is given by ([I]). The path integral is to include all paths satisfying rii{($) = rij(O) 
and 4>i(0) = 4>i{(3). The variables are integers and therefore the path integral has to be 
performed stepwize [^TTj. We limit ourselves to zero temperature, i.e., (3 — oo. 

The first step of the derivation follows closely previous works pOf , and is therefore given 



in Appendix A. Using the Villain approximation and the duality transformation method, 
the path integral over the charge and phase degrees of freedom, and is transformed 
to a path integral over a 3-component integer field J™ r describing the vortex charge and 
density, and a 3-component real gauge field /Q, to which J™ r is coupled. This gauge field 
describes the charge degrees of freedom, to which it is related through its derivatives. The 
field strengths associated with this gauge field, ^r&^d^Ki^ are the Cooper-pairs current 
and density on the i'th site. In terms of J vor and /C, and in a gauge in which A • /C = the 
action is given by, 
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S vor = J2i {i(pT - + iJr ■ (£ + R xt ) + s^Wjii^) 2 + (AM 2 ]} 

(3) 

+^E iJ (A x /Q^A x /C,) 

where A x }C ext = 2nhn x . This action describes the vortices as bosons on a lattice, interacting 
with an externally put gauge field K ext whose spatial curl is a constant, given by 2irhn x , 
and with a dynamical gauge field K. As expected from the similarity between Magnus and 
Lorenz forces, a moving vortex is affected by the Cooper-pairs on the dots in the same 
manner a charged particle is affected by a magnetic field. Moreover, the Josephson currents 
between the dots affect the vortices in the same way an electric field affects charged particles. 
The last three terms of the action include the self energy of the field /C. They are simply 
understood once the relation between K, and the Cooper-pairs currents and densities is 
taken into account. The first two are the kinetic energy of the Josephson currents (the 
transverse part of that current is ^A/Co, and -%/C is the longitudinal part). The last term 
is the charging energy (the net charge on the z'th dot is A x /Q). The transverse part 
of the current satisfies a two dimensional Gauss law A 2 /Co = A-K 2 Ejp vor and mediates a 
logarithmic interaction between the vortices. The excitation spectrum of JC is the spectrum 
of longitudinal oscillations of the Cooper-pairs, i.e., the plasma spectrum of the array. 

Our next step is a formulation of the continuum limit of the action @. When doing that, 
two points should be handled carefully. The first is the periodic potential exerted by the 
lattice on the vortex. This potential was studied in detail by Lobb, Abraham and Tinkham 



22]. Since Currents do not flow uniformly within the array, the energy cost associated with 



a creation of a vortex depends on the position of its center within a plaquette (i.e., on the 
precise distribution of the currents circulating its core). This energy cost is periodic with 
respect to a lattice spacing of the array, and is independent of the sign of the vorticity. 
The origin of this potential can be visualized using the analogy with 2D electrostatics. In 
that analogy, a vortex is analogous to a charge in a two-dimensional world. A vortex on 
a lattice is then analogous to a charge in a two dimensional world in which the dielectric 
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constant varies periodically with position. The electrostatic energy of such a charge varies 
periodically with position, too, and is independent of the sign of the charge. This energy 
cost can then be interpreted as a periodic potential exerted by the lattice on the vortices. 
The characteristic energy scale for that potential is Ej. Its amplitude and functional form 
were studied in Ref. P^fl . The amplitude was found to be 0.2Ej and 0.05Ej for square and 
triangular lattices, respectively. 

A convenient way to incorporate the periodic dependence of the vortex potential energy 
on the position of the vortex core within a plaquette is by replacing the Josephson energy 
Ej in the action @ by a periodically space dependent function ej(f), that is non-zero only 
along lattice bonds. The period of ej(r) is obviously the lattice spacing. The energy cost 
involved with the Josephson currents then becomes / d f Sn2 ^^ [(V/C (r, t)) 2 + (JC(r, t)) 2 ], 
and that energy cost confines the currents to the lattice bonds. The effect of the spatial 
dependence of ej(r) on the interactions mediated by /Co,/C is discussed below. 

The second point to be handled carefully when transforming to a a continuum action is 
the short distance cut-off on the capacitance matrix. The model we employ does not attempt 
to describe statics and dynamics of Cooper-pairs within superconducting dots. Thus, its 
continuum version should not allow for excitations of K at wavelengths shorter than the 
lattice spaing. This constraint is taken into account by introducing a high wave-vector 
cut-off to the capacitance matrix, as it was done in Ref. [14]. 

Taking into account the two points discussed above, the continuum limit of the action 

(I) is, 

S =JdtJdf\ i[p vor (f, t) - n v ]JC {f, t) + if(f, t) ■ [/C(f, t) + IC ext (r)} 



^[(VK (r,t)) 2 + (t(f,t)) 2 }\ (4) 
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+^fdtfdrfdf* [V x IC(r,t)\ C-\r-T f ) [V x 1C(f \t) 
When the fields /Co,/C are integrated out they mediate mutual interactions between 
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vortices and self interactions of a vortex with itself. The spatial dependence of ej(r) does 
not significantly affect mutual interactions between vortices whose distance is much larger 
than one lattice spacing. It does, however, potentially affect the self interaction. 

Consider the action associated with a single vortex. As discussed above, due to the 
spatial dependence of ej(r), the interaction of the vortex with /C yields a periodically space 
dependent potential energy f22|. The coupling to /C results in a kinetic energy. To see that, 



note that the vortex current corresponding to a moving vortex whose center is at fo(t) is 
jvor _ f Q (j-^fi^_ r (t)). Substituting this expression in the action (P, we find the part of 
the action that depends on the vortex velocity, r , and the gauge field it interacts with, /C, 
to be 

iJdt? (t)-)C(r ,t) 

+ fdtU df^j^(t(f, t)f + dfj dr 1 [V x /C(f, t)] C-\r- f) [V X JC(f>, t)} j 

(5) 

The gauge field /C can be integrated out, with the resulting effective action for the vortex 



velocity r being non-local in time. However, as pointed out by Eckern and Schmid [14 



and by Fazio et.al. ||20|| , the time non-locality can be neglected as long as the characteristic 
frequencies involved in r (t) are smaller than the Josephson plasma frequency ^\/8EjEc- 
This neglect is possible due to the gap in the excitation spectrum of /C, a gap that makes 



the time non-locality short ranged ||23|| . Having neglected the time non-locality, we find the 



effective action for the vortex velocity Tq to be, 



1 • 2 

dt-m bare f (t) (6) 



where rribare, the vortex bare mass, is defined by rribare = tw~- [0] Thus, the interaction of 



4E C 

the vortex with K, results in a kinetic term. 

Having integrated out both /Co and /C we have turned the single vortex action into 
an action of a charged particle interacting with a periodic potential and a magnetic field 
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2nhn x . For ^ < 1 the effect of the periodic potential is to change the bare mass into 
an effective band mass. Since the lowest band is the relevant one for bosons, the band 
mass is always larger than the bare one p4|. The effective band mass for n x = was 



studied both theoretically and numerically by Geigenmuller [^5J and by Fazio et.al. pO 
(see also references therein). While a qualitative estimate of the mass is easy to arrive at, 
a quantitative determination depends on the precise details of the periodic potential, and 
is therefore hard to obtain. Qualitatively, the tight binding limit, in which Ej ^> Ec, is 
distinguished from the weak periodic potential limit, in which the opposite condition applies. 
In the former, the effective band mass is 



m band ~ ^J^e^ (7) 



where are numbers of order unity p0[ [25|. In the latter, 

E 



m band ~ m bare ( 1 + {^-f ) (8) 



The regime of parameters we are interested in lies between the two limits. It is therefore 
reasonable to assume that the band mass is larger than, but of the order of, the bare one. 

We now turn to discuss the many vortices configuration. As we have argued above, the 
discreteness of the array does not significantly affect the mutual interaction between vortices. 
Therefore, for the study of this interaction we may replace ej(r) by Ej. Then, the action 
(f|) can be written in momentum space as 

S = fdtJjgp fa- ? ~ n v 5{q)]K^+iJ}° r ■ + s^jtl^Cod 2 + l4 



1 21 



(9) 

where p^ or , J| or , K.^, /C^ 1 '*, C~ l (q) are the Fourier transforms of the corresponding quantities. 
The momentum space representation used in Eq. (|9]) is convenient for the integration of the 
field /C. The integration out of the time component, /Co, yields a density-density interaction 
between the vortices, of the form | J dq^\p^ or — n v 8(q)\ 2 , where q = \q\. When transformed 
back to real space, this interaction is 

13 



ttEjJ dfj dr"(p vor (r) -n v )\og\f-r'\(p vor (f*) -n v ) (10) 

Similarly, the integration of /C yields a current-current interaction between the vortices. In 
the gauge we use, K,^ has only transverse components. It therefore mediates an interaction 
only between the transverse component of the vortex currents. Integrating out /C, and 
neglecting again the slight non-locality in time, we find that the current-current interaction 
is given, in momentum space, by 

Tr ' <U1 'JZ? = ttet I "dgWI 2 ( n ) 



where J v f^ = - x JH or is the transverse component of J| or , and the integral over q is cut-off 
at q — 2tt. The current-current interaction is described in real space as, 



n 2 



fdf f dr> ( dq[J vor {r) x q\ [ J vor (j*) x q\— e ^~^ (12) 
16Ec J J J q 2 



As pointed out in Ref. |L4]], this current-current interaction includes both a self interaction 
term, that assigns a mass mb are to each vortex, and a velocity-velocity interaction between 
different vortices. The former was discussed in the context of a single vortex. The latter is 
short ranged. For large separations, it is inversly proportional to the square of the distance 
between the interacting vortices. 

Eqs. ( |i~0"|) and (|1^) both neglected the effect of the lattice structure of the array on the 
vortex-vortex interactions. Similar to the common practice in the analysis of electrons on a 
lattice, we assume that the sole effect of the lattice is to modify the single vortex mass from 
the bare mass rribare to the effective band mass rriband- The current -current interaction in 
Eq. (|12D includes a self interaction that assigns a mass vribare to each vortex. To account for 
the modification of the mass by the lattice, we add another kinetic term to the action, of 
the form M* J d f^ pV0 ^L , where M* = rriband — rnbare- Altogether, then, the effective vortices 
action becomes, 
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S wr (J wr ) = fdt\ J df\M*^^ + tJ vor {r) ■ K, ex \r) 



fdrfd? Jdq\J vor (r) x q\[J vor (r f ) x q] 



(13) 



+nEj(p vor (r) -ri v )log\r-r'\(p vor (r') -n v ) 



Eq. (0) is a concise description of the dynamics of the vortices, since the only dinamical 
fields it includes are those of the vortices. It describes the vortices as interacting particles 
of mass rriband and an average density n v , under the effect of a "magnetic field" 2irhn x . The 
vortices "filling factor" is then indeed — . 

& n x 

The current -current interaction term in the action (|T3|) is somewhat inconvenient for 
calculations. Thus, in our analysis of the quantum Hall fluid of vortices in the next section 
we choose to reintroduce /C, and consider the vortices as particles of mass M* interacting 
with a dynamical vector potential /C, as well as with K ext and with one another. 

We conclude this section with a few remarks regarding the dependence of its results on 
the form of the capacitance matrix. The capacitance matrix determines the bare mass of the 
vortex (see Eqs. (|5p and (Q)) and the form of the vortex current-current interaction (see Eq. 
(|l2l)). So far we considered a capacitance matrix that includes only nearest neighbors cou- 
pling. The inverse capacitance matrix describes then a two dimensional Coulomb interaction 
between Cooper-pairs on the superconducting dots. In Fourier space, it is proportional to 
Inclusion of capacitances to the ground and/or capacitance between dots that are not 
nearest neighbors result in a screening of that interaction. Then, at small q, C~ l (q) oc q~ a , 
with < a < 2. This screening has two consequences. First, the kinetic energy cost in- 
volved in a vortex motion, i.e., its bare mass, is affected. Second, the excitation spectrum 
of K, becomes gapless. We now examine these consequences. Consider a vortex moving in a 
constant velocity v. As seen from Eq. (|3p, a moving vortex acts like a source for the vector 
potential /C. The (imaginary time) wave equation for JC can be derived from Eq. (|5|). In 
Fourier space its solution is, 
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where v± = is the transverse part of the velocity vector. The longitudinal component 
of JC vanishes in the gauge we use. The kinetic energy cost associated with the motion of a 
vortex is the energy cost of the fields K. and V x K, it creates. It is composed of two parts. 
The first, / d f^l Ej K , is the kinetic energy cost of the longitudinal currents created by the 
motion of the vortex. The second, e 2 2 3 J drj dr 1 V x K{f) C~ l (r — ?) Vx /C(r*) is the 
cost in charging energy. A moving vortex induces voltage drops between the superconducting 
dots, and those result in a charging energy cost, determined by the capacitance matrix. The 
first energy cost is proportional to i> 4 , while the second is proportional to v 2 . Thus, the 
bare mass is determined by the charging energy. Transforming Eq. ([TJ]) to real space, and 
substituting into the expression for the charging energy, we observe that the charging energy 
is finite as long as a > 0, and diverges logarithmically with the system size for a = 0. The 
gapless excitations of /C, characteristic of a < 2, play a role when vortices accelarate or 
decelarate. Then, the coupling of the vortices to these excitations (the "spin waves" p0| ) 
becomes a weak mechanism for dissipation of a vortex kinetic energy [TJj]. For the present 
context we note that the effect of a weak dissipative mechanism on the quantized Hall effect 



was studied by Hanna and Lee [26]. While some properties of the effect are affected by such 



a mechanism, its main features are not. 



4. The quantum Hall fluid of vortices 

4.1 General discussion 

In the previous section we established the mapping of the vortex dynamics in a Josephson 
junction array on the problem of interacting charged particles in a magnetic field. We have 
also identified the ratio — as the vortices filling factor. In this section we examine the 

Tlx & 

formation of QHE fluid state of vortices at appropriate filling factors. While so far we have 
emphasized the similarities between the dynamics of the vortices and that of electrons in 
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a magnetic field, in this section we must study the differences between the two. We begin 
by a general discussion of two of the differences. Then, we turn in the next subsection 
to a detailed calculation, using the Chern Simon Landau Ginzburg approach to the QHE, 
developed by Zhang, Hansson, Kivelson and Lee [^7] . 



The first difference is in the statistics: while vortices are bosons, electrons are fermions. 
This difference changes the values of the "magic" filling factors, and eliminates the possibility 
of a QHE in the absence of interactions. The filling factors at which bosons form quantum 
Hall fluids are -, where p, q are integers, and one of them is even [28fl. Fermi liquids of the 
type discussed by Halperin, Lee and Read form at filling factors r 2 n+i) ' where n is an 
integer. 

The second difference is in the interaction: the logarithmic interaction between vortices 
is of longer range than the Coulomb interaction between electrons. This difference leads to 
a modification of the quantized Hall conductance, a modification of the charge of Laughlin's 
quasiparticles, and, perhaps most interestingly, to a modification of one of the diagonal 
elements of the linear response function. While for a short range interaction these diagonal 
elements vanish in the long wavelength low frequency limit (q, to — > 0), reflecting the lack of 
longitudinal dc conductance in the QHE state, we find that the logarithmic interaction makes 
one of the diagonal elements non-zero. In fact, rather than describing insulator-type zero 
longitudinal response, as expected from a QHE system, this element describes a longitudinal 
response of the type usually associated with a superconductor. These consequences of the 
logarithmic interaction are all derived in detail in the next subsection, where we also study 
the difference between the linear response function and the conductivity. In this subsection 
we preceed the derivation by a discussion of a thought experiment that makes the role of 
the logarithmic interaction physically transparent. The thought experiment we consider 



was extensively used in the study of the Quantized Hall Effect, e.g., by Laughlin |30[ and 
Halperin ||31|| , and was proved very useful in understanding various aspects of the effect. 

Consider a "conventional" electronic quantum Hall system, in which a disk shaped two 
dimensional electron gas (2DEG) is put in a strong magnetic field, and a thin solenoid 
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threads the disk at its center. The flux through the solenoid is time dependent, and is 
denoted by $(£). If the electrons on the disk are in a QHE state, the current density at any 
point is perpendicular to the total electric field at that point. The time dependence of the 
flux induces an electric field in the azymuthal direction, given by where c is the 

speed of light, and r is the distance from the center. Due to the finite Hall conductance, 
this electric field induces a radial current, and, consequently, a charge accumulation at the 
center of the disk. The charge accumulated at the center during the interval < t < t is 
given by ^§^-($(to) — $(0)) (where a xy is the dimensionless Hall conductivity and <3> = is 
the flux quantum). This charge accumulation, in turn, creates a radial electric field. Now, 
if the electrons interact via a Coulomb interaction, the radial electric field is proportional to 
■\, i.e., it decays faster than the azymuthal one. Then, far away from the center the electric 
field is predominantly azimuthal, and the currents are predominantly radial. However, if the 
electrons interact via a logarithmic interaction, both the radial and azymuthal components 
of the electric field are inversly proportional to r, and thus their ratio is independent of r. 
The current then has both radial and azimuthal components, and their ratio is independent 
of r, too. Moreover, the azymuthal component of the current is proportional to the flux at 
the center, and not to its time derivative. 

The two components of the current and the charge accumulated in the center can be 
calculated using classical equations of motion, since in the absence of impurities, the classical 
and quantum mechanical calculations coincide. Consider, therefore, the hydrodynamical 
equation of motion of a fluid of electrons in a magnetic field, whose electronic density and 
velocity fields are denoted by p(f) and v(f), respectively. Assuming a uniform positive 
background charge density p on the disk, this equation of motion is 

mp(r)i?(r) = -p{r)v{r) x B- J d?VV e _ e {r - T*)(p{f) - p)p(r) + ^Pif) ( 15 ) 

where v(f) is the complete time derivative of the velocity field, V e - e is the electron-electron 
interaction potential, and we use a system of units where e = c = 1. The initial conditions 
corresponding to the scenario discussed in the previous paragraph are p(r) = p, v(f) = 
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and <&(t = 0) = 0. Due to the circular symmetry of both Eq. flT5|) and its initial and 
boundary conditions, the current and density remain circularly symmetric when the flux 
is turned on, and the electron-electron interaction term can be written as V^_ e (r)Q(r) p(r) 
where Q(r) = Jq dr'2-Kr'(p(r') — ~p) is the net charge within a distance r from the origin, and 
V r e '_ e (r) = d^"^- 1 . The conservation of charge constraint implies Q(r) = —2nrp(r)v r (r), 
where v r is the radial component of v. The azymuthal component of Eq. (|15|) can therefore 
be written as, 

2irrp(r)v4 = - u c Q (16) 

m 

with uj c = For values of r far away from the center but not close to the edge the 
density p(r) remains approximately equal to p all along the process, and thus Q{r) is r- 
independent. Within that approximation, and for such values of r, the azymuthal equation 
can be integrated and substituted in the radial one. The latter then becomes, 

u 2 c Q p uj c <& 1 , 1_ 2 

/Wr = o 7^ V e-e( r )Q( r )P + -p v A r ) ( 17 ) 

27rr m 27rr m r 

where the last term is the centrifugal force. Suppose now that the flux $ is turned adiabat- 
ically on from zero to $(to) i n the interval < t < t . For times t ^> to the velocity field is 
purely azimuthal and v r = 0. Then, if the potential gradient W e - e decays faster than -, so 
does also the azymuthal velocity v^, and 

= ^ (18) 

where v = -|$o is the filling factor. If $(to) — then the charge accumulated near the 
origin is the charge of Laughlin's quasiparticle, namely —ev. 

However, in the case of a logarithmic interaction, V r e _ e (r) = — Vologr, 



(19) 



Q 



hujc 
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The azymuthal current is then indeed inversly proportional to the distance from the origin, 
and proportional to the flux $. In fact, the current is related to the vector potential created 
by the solenoid, A so1 , via a London-type equation 

V x J — v x A so1 (20) 

m hu c + Vqv 

Thus, the longitudinal response of the electrons on the disk to the vector potential created 
by the solenoid resembles the longitudinal response of a two dimensional superconducting 
disk in a similar situation. 

Two conclusions can be drawn from the above thought experiment. First, for logarith- 
mically interacting particles, the charge of the Laughlin quasiparticle does not equal the 
filling factor, but depends on the interaction. And second, the transverse part of the linear 
response function, relating a transverse current to an externally applied transverse vector 
potential, resembles that of a superconductor. Since this response function is proportional 
to the transverse current-current correlation function, the latter should be expected to re- 
semble a supeconductor, too. Both conclusions are substantiated in the next subsection, 
and are applied to the study of the quantum Hall fluid of vortices. 



4.2 A study of the vortices QHE state by the Chern Simon Landau 
Ginzburg approach 

In this subsection we use the Chern Simons Landau Ginzburg approach to further analyze 
the properties of the quantized Hall fluid of vortices formed at appropriate values of — . We 
start by writing a Landau-Ginzburg action that describes the dynamics of the vortices. We 
then perform a Chern-Simons singular gauge transformation that attaches an even number 
of fictitious Cooper pairs to each vortex. The resulting action, in which the order parameter 
describes transformed "composite" bosons, is convenient for a saddle point analysis. We 
find the uniform density saddle point that describes a superfluid of composite bosons. By 
expanding the action to quadratic order around that saddle point we calculate the response 
function of the vortices to an external probing field. This response function, denoted by E, 
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is the ratio of the vortex density and current J vor to an infinitesimal probing field /C p applied 
externally to the system. The matrix S is calculated by adding an external probing field /C p 
to the Lagrangian, and integrating out all the other fields to obtain an effective Lagrangian 
L e ff in terms of /C p only The three components of J vor are then given by 
and the elements of S are W> 



dL e ff 



K*=0 



_d d_ 



L e// (/C p ) 



(21) 



/cp=o 



The physical meaning of S, as well as the important distinction between the response to /C p 
and the response to the total field /C p + /C, are discussed after the calculation is presented. 

The Landau-Ginzburg action that describes the properties of the vortices as they were 
found in section (3) is, 

Slg$,1C) = Jdt\ fdr[nrd t 4> + ^\(ihV - K** - K)$\ 2 + 



+ / dfj df TrEj(\^(f)\ 2 — n v ) log |r — r l \(\ip{r')\ 2 — n v ) 



(22) 



The fields i/j,K. ext ,)C all depend on r and on t. For the brevity of the expressions we omit 
this dependence whenever this omission does not lead to confusion. The field the order 
parameter for the vortices, satisfies bosonic commutation relations. 

Restricting our derivation to the "fundamental" fractions -, where rj is an even number, 
our first step in analysing the action (p2|) is the Chern Simon singular gauge transformation, 
in which the field i/j(r,t) is transformed to 



(23) 



where arg(r — f) is the angle the vector f — r 1 forms with the x axis. Since rj is an even 
integer the field if) has the same statistics as ■0, i.e., bosonic. Note that \ip(f, t)\ = \tp(f,t)\. 
The singular gauge transformation shifts the phase of the field. Denoting the phase of 
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by 8(6), the Chern-Simons transformation amounts to V8(r,t) = V8(r,t) — ~a(r,t) where 
a(r,t) = fir] J df^^y^^-\ijj(r , ,t)\ 2 . The Chern-Simons field a has a gauge freedom, which 



we fix below. Writing ip(f,t) = \Jn v (r,t)e %e ^ r,t \ the above Landau-Ginzburg functional 
becomes 

S LG (n v , 6, K, a) = / dtl J dr\in v (hd t 8 - a ) + ^(hVd - lC ext -K + a) 2 + 



(24) 



+ / drj dr 1 nEj(n v (r) — n v ) log \ f — r'^n^r' 



+^[V x £(r)}C-\r- r-)[V x /C(r')] 

While the Euler-Lagrange equations of motion for the original field ip required its phase to 
be multiply valued (for a minimization of the kinetic energy), the equations of motion for 
the transformed field ip allow for a solution in which the magntiude and the phase of the 
field are constant. It is straight forward to see that the action ( p4[) is minimized when, 



n v (f,t) = n v 



(25) 



V0(f, t) = 



a(f, t) = IC ext (r) 



(26) 



(27) 



K(r, t) = (28) 

This minimum of the action describes a state in which the vortex density is constant on 
the average and the Chern-Simons field a cancels K ext on the average. We now expand 
the action around these minimum values. Around the saddle point the phase 8 is singly 
valued, and thus we can choose a gauge in which 8 = identically, and the field ip is real. 
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Writing, then, n v = n v + 5n v , ip = Jn^ + and a = IC ext + 5a, we find that the quadratic 



deviations from the extermum point fl2"5f )-(p8|) are described by the following action, 



S LG (5n v ,JC,a) « / dt / dr[ia 5n v + + 5a) 2 + ^ 



•2 



+ fdtfdrfdr , \'irEj6n v (r) \og\r - f*\5n v (f*) + ^(V x /C^C^f - r*)(V x 



(29) 



For the calculation of the conductivity matrix and the correlations functions, we imagine 
coupling the vortices to an infinitesimal 3- vector probing field /C p . Naturally, linear response 
functions are more conveniently described in Fourier space. Thus, we write the action (|29|) 
in the presence of the probing field /C p , in Fourier space, as 



(2tt)2 



i(oo(q) + /Cg(q))5n,(-q) + fe\KP(q) + /C(q) + 5a(q 



h 2 \qSn v (q) 

" 8M* n„ 



+ 8^«q)| 2 + ^ v °a^)q v aM) 



f ^|^(q)| 2 + x /^(q)! 2 ^ 1 ^ 



(30) 



where qo = a>. We choose the Coulomb gauge for ZC P , i.e., q ■ )C p (q) = 0. Thus, the probing 
field becomes a 2-component vector. Since the three components of J vor are constrained by 
the conservation of vorticity, J vor is effectively a two-component vetor, too, and S a/ g is a 
2x2 matrix. The indices a and f3 take the values (for the time component) and _L (for 
the component perpendicular to q). The integration of the fields 5n v ,K.,a is easily carried 
out, since (|30|) is quadratic in all fields. The resulting effective Lagrangian is, 



(K».) 2 



L e// (K?) 



M*u 2 _|_ 2-kEj _|_ hq A 
2nhr]n v hr] 8nM*n v r} 



+ 



2D 



iqJC p JC p ± 



(31) 



where D = 27r7w7 



Ml _i_ 



D = 2?T^ 



Ai* 



+ 2m, 



fcare 



1 



In the limit of q, u> — > and for (7 (g) = C nn q 2 
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Eq. (|2T| ) expresses the matrix elements of E in terms of second derivatives of this effective 
Lagrangian with respect to /Cq and /C^. Each of the four components of the matrix E Q/ g 
warrants a short discussion. First, the transverse component of the vortex current and the 
transverse component of the gauge field /C p are related, in the limit q, lo — > 0, by 

= ZkEjD 

1 An 2 hr]Ej + 2nh 2 rfD 1 K ' 

This London-type of relation was anticipated by Eq. (pO]). It is characteristic of supercon- 
ductors, and is very different from the insulating behaviour characteristic of the diagonal 
components of the response functions in QHE systems. This difference results from the 
static vortex-vortex interaction being of a long range. Defering the discussion of the effect 
of Eq. (|32|) on the longitudinal conductivity to a later stage, we now point out its effect on 
the vortex current-current correlation function. By the fluctuation-dissipation theorem, 

(JT/r) ?> = Im£ ±± (?>) 

(33) 

/do/P(^)(Jrj r ) 9> , = Re£ ±± (g» 

where P denotes the principal part of the integral, and the second line is an application 
of Kramers-Kronig relations f32| . For an insulator, Re£((f, to) oc lo 2 when lo — > 0. This 



is also the case for QHE systems with short range interactions |33 |. For a superconductor 



ReS(g, lo) approaches a constant in the lo — > limit. As we now see, so is also the case for a 
QHE system in which the interactions are logarithmic. In the particular problem we study, 
this constant is — , 2t It-tvts 2 n - 

Second, we note that the compressibility of the vortex fluid vanishes in the limit q, lo —>■ 0, 
as is manifested by the absence of low frequency poles in the density-density correlation 
function. Like its electronic analog, the quantum Hall fluid of vortices is incompressible. 

Third, the Hall component of the linear response function is given, in the limit q, lo — >• 0, 

by, 

S °' ± = ^ + 2nh V (34) 
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If the vortex-vortex interaction was of shorter range, the long wavelength limit of Eq,± would 



satisfy S c 



2-Khrj 



seemingly demonstrating the quantization of the Hall conductivity |]27 



We further comment on the difference between the two expressions below. 

The qualitative effect the logarithmic interaction has on the transverse and Hall com- 
ponents of the linear response function raises the following question: does the conductivity 
of the system we study have the properties of the conductivity matrix of a QHE system, 
namely, zero longitudinal conductivity and quantized Hall conductivity? To answer this 
question, we clarify the relation between the response function E a/ g and the vortex conduc- 
tivity matrix. A similar relation was discussed, in the context of the QHE, by Halperin 
|3lfl, Laughlin [34j], Halperin Lee and Read |29| and Simon and Halperin [35|| . The transport 
of vortices in the array is probed by externally applied (number) density and current of 
Cooper-pairs, given, respectively, by x /C p , and — ^ V/Cq — 2^K^ P - ^he ma trix £ is 

defined such that J vor = S/C p . However, the vortices themselves contribute to the Cooper- 
pair density and current, with the most trivial contribution being the circulation of current 
around each vortex center. The total Cooper-pair density and current are therefore given by 
the derivatives of a total gauge field, composed of the probing field /C p and the field induced 
by the vortex density and current, denoted by /C md . The latter is proportional to the vortex 
density and current J vor . Thus, we can define a matrix V such that 

/C ind = V3 vor = V±rC p 



Consequently, the total field is /C tot = (1 + Vt)KP, and 

j vor = tii + vty 1 ^ 

Thus, the matrix S(l + VT,)^ 1 relates /C tot to the vector 



p uu, 

jvor 
\J± J 



(35) 
(36) 

The vortex conduc- 



tivity matrix a vor , relating the vortex current 



TV 
J \\ 



to the total driving force vector 



■ q y-tot 

, ~i^rC tot 
\ 2-kH 



is then, 
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q 





t(l + Vt) 



-i 



/ z 2^ Q \ 

q 







(37) 



/ 



where the leftmost matrix converts p vor to J|'| 0r . 

Eq. fl37|) defines the vortex conductivity matrix in terms of the matrices S and V\ The 
matrix £ is defined by Eqs. ( pip and (|3T]) . The matrix V, relating J vor to JC ind , is specified 
by the action (EDI) to be 



V 



47T 2 gj 











(38) 



% + -^q-2C'- 1 (q) ] 



The upper left element describes the field /Co created by a vortex density p vor . The gradient 
of that field is the transverse current circulating around the vortex center. The bottom right 
element describes the field JC± created by a transverse vortex current J^ or , and is obtained 
from Eq. ( |3~0l) by taking its derivative with respect to /Cj_(q). 

Substituting the matrices S and V to Eq. ([37]), we find, to leading order in q,u, 



I 



; mM* 1 
' r) 2 n v r) 



\ 



_1 

V n 



•ojM* 



(39) 



ri 2 n v ) 

In the limit q, to — > the diagonal terms vanish, and the vortex current satisfies 



J v 



z x (qlCo + uK) 



(40) 



Eqs. (^) and ([!(]) describe a quantized Hall effect: the current is purely perpendicular to 
the total "driving force", and the Hall conductivity is quantized. Contrasting Eqs. (|3^) and 
(|34]) with Eq. (f40"D we can finally summarize the effect of the logarithmic interaction on the 
linear response of the system: the dc conductivity, which is the q, uj — > response to the total 
driving force, has the usual form of the quantum Hall conductivity, and is unaffected by the 
interaction. The correlation functions, on the other hand, determined by the response to the 
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externally applied driving force, are affected by the interactions even in the q,uj — > limit, 
with the most notable effect being on the transverse current-current correlation function. 

We conclude this section by relating the vortices conductivity, calculated above, to the 
electric conductivity and resistivity, which are the quantities typically measured in experi- 
ments. The electric conductivity is the matrix relating voltage drops (or, in the continuum 
limit, electric fields) between superconducting dots to the electric Josephson current flowing 
in the array. The electrostatic potential at a point r is given by -% J dr*C~ 1 (r — r')'V x JC(f f ). 
Thus, in Fourier space the q component of the electrostatic potential is i^qC^ 1 (—q)lC±(q) 
and the longitudinal electric field is —^q 2 C~ 1 (—q)IC±(q). Now, by deriving an equation of 
motion for )C± from the action we see that a dc transverse vortex current J v f r creates a 
field /Cj. given by 

J% = -^fC- 1 ®^® (41) 
i.e., a transverse vortex current J v °~ induces a longitudinal electric field ^Jj°^ ■ A similar 
argument regarding the relation of the longitudinal vortex current to the transverse electric 
field leads to the conclusion that a vortex current JH or creates an electric field — z x Jy°J. The 
Josephson charge current, on the other hand, is —i^z x (qTCo + Lo)C±), i.e., it is proportional 
to the "driving force" acting on the vortices. Thus, the matrix relating the Josephson 
current to the electric field is proportional to the matrix relating the driving force acting on 
the vortices to the vortex current, or, explicitly, 

f-w*- (42) 

where p el is the electric resistivity matrix of the array @. This result can be simply concluded 
from Eq. (0). The right hand side of that equation is the Josephson current, divided 
by 2e. The left hand side is proportional and perpendicular to the electric field. The 
electric field is then proportional and perpendicular to the Josephson current, with the 
proportionality constant being tJ^w^- The quantum Hall fluid of vortices manifests itself in 
electronic properties of the array - the longitudinal electric resistivity vanishes, and the Hall 
electric resistivity is quantized. 
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5. Conclusions 

In the previous sections we presented a study of the transport of vortices in an array 
of Josephson junctions described by the Hamiltonian (|1|). In particular, we focused on a 
quantum Hall fluid formed by the vortices at appropriate values of In this section we 
summarize the results of this study, and comment on a few open questions. 

Our study was motivated by the analogy between Magnus force acting on vortices and 
Lorenz force acting on charges in a magnetic field. In this analogy, fluid density plays a 
role analogous to a magnetic field, and fluid current density plays a role analogous to an 
electric field. Quantum mechanics extends the analogy further: a fluid particle is found to 
play a role analogous to that of a flux quantum. This analogy motivates the search for a 
quantized Hall effect for the vortices. The vortices' filling factor is identified with the ratio 
of the vortex density to the fluid density. This ratio is very small for superconducting films, 
and it is this smallness that motivates the study of the Josephson junction array. Due to 
the periodicty of the spectrum of the Hamiltonian ([]]) with respect to the parameter n x , the 
effective filling factor becomes - ^ odl ^ , which can be made of order unity. 

The dynamics of the vortices in a Josephson junction array was studied in section (3). It 
is found to be that of massive interacting charged particles under the effect of a magnetic field 
and a periodic potential. The magnetic field is 2nhn x . The effect of the periodic potential is 
taken into account in an effective mass approximation, changing the mass from a bare mass 
to an effective band mass. The effective mass is exponentially large for Ej 3> Ec, and of 
the order of for Ej>Ec- Being interested in a phenomenon resulting from a motion of 
vortices, we obviously consider the latter regime. The mutual interaction between vortices 
consists of a velocity independent logarithmic interaction, whose strength is proportional to 
Ej, and a short ranged velocity- velocity interaction. 

In view of the mapping of the dynamics of the vortices on that of massive interacting 
charged particles in a magnetic field, the existence of a quantum Hall fluid phase is to 
be expected. In Section (4) we examine some of the properties of that phase, but we leave 
unexplored some other important propoerties. Most notable among the latter are the regime 
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of vortices filling factors at which the quantum Hall fluid is the lowest energy state, and the 
energy gap for excitations above that fluid. 

Our study of the quantum Hall fluid is performed by means of the Chern Simon Landau 
Ginzburg approach to the quantum Hall effect. When ^ = - with 77 being an even integer, 
the vortices Landau-Ginzburg action is found to have a saddle point corresponding to a 
quantum Hall fluid. By hierarchical construction such saddle points can be found for ^ = |, 
with p, q being one even and one odd integer. The properties of the vortices quantum Hall 
fluid are studied within a quadratic expansion of the action around the corresponding saddle 
point. We find that the vortices conductivity matrix shows a typical QHE behavior, i.e., zero 
diagonal elements and quantized non-diagonal elements. However, we find the q, uj — > limit 
of the current-current correlation functions in the ground state to be different from those of 
a typical quantum Hall state, due to the long range logarithmic interaction. In particular, 
the transverse current-current correlation function is predicted to bahave like that of a 
superconductor, rather than an insulator. For large arrays (larger than an effective London 
penetration length) the vortex-vortex interaction is screened. Then, both the conductivity 
matrix and the correlation functions are expected to behave, in the dc limit, like those of a 
typical quantum Hall state. 

A necessary condition for the quantum Hall fluid to be the lowest energy state is pre- 
sumably that the ground state at n x = and n v 7^ (infinite filling factor) is a superfluid 
of vortices, i.e., an insulator. The observation, by van der Zant et. al. ||15|| , of a magnetic 
field tuned transition points at the regime of parameters in which this condition is satisfied, 
namely, Ej rs Eq and 0.3 > n v > 0.15. In this regime of parameters we expect the quantum 
Hall fluid to be the ground state at large filling factors, and the Abrikosov lattice to be the 
ground state at small filling factors. This expectation is based on the phase diagram of a 
two dimensional electron gas. For the latter, if the ground state at zero magnetic field is a 
Fermi liquid, then the ground state at large filling factors (>0.2) is the quantum Hall fluid 
and the ground state at low filling factors is the Wigner lattice. 
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Appendix A — The Villain approximation and the duality transfor- 
mation 

The starting point of this appendix is the expression of the partition function as a path 
integral over the phase and number sets of variables {<pi}, {fii}, Eq. ([|). Using the Villain 
approximation and the duality transformation we transform that path integral to a path 
integral over an integer 3-component vector field J vor , describing the vortex density and 
current, and a real 3-component vector gauge field, /C. The action in terms of J vor and /C, 
to be derived below, is given by Eq. (H). The following derivation follows the method of 
Fazio et.al. |2(J 



In the Villain approximation the imaginary time integral is done in discrete steps, where 
the size of each step, denoted by r , is of the order of the inverse Josephson plasma fre- 
quency ujj = %~ 1 \/9>EjEc- Each term in the Josephson energy part of the path integral 
is approximated by a Villain form (we put % = 1 throughout the appendix, and restore its 
value in the final formula), 

e — r Ej(l-cos(^— Aij)) y^oo g _ I TQ £; J (0. J ._ Aij+2-KVij) 2 



1-1 



,. (43) 

-c^oo / rp r j -^mT+ l PijTo{<t>ij~A l: j+2iTv i: j) 



This approximation is valid for EjTq>1, and gets better as EjTq gets larger. However, it 
retains the most important feature of the Josephson energy, the periodicity with respect to 
(j), for all values of Ejt . Alogether, then, the approximation we discuss holds for Ej > Eq. 
The significance of the field V{ can be understood by noting that A x Vi describes the density 



of vortices |36| . As for the real variable pij, as shown below, it describes the Josephson 
current along the bond ij. Since the Josephson energy includes a sum over all lattice bonds, 
the Villain approximation introduces a variable pij to each lattice bond. Thus, we can regard 
p as a vector defined for each lattice site, such that p ix corresponds to the bond i, x and p iy 
corresponds to the bond i,y. Similarly, the difference <pi — <pj, the integral J? A ■ dl and the 
variables can be regarded as vectors A^j, Ai and V{. 

Next we apply the Poisson resummation formula to the dependent part of the action. 
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By doing that we make the n« variables real numbers rather than integers and add a time 
component to the integer- valued vector field Uj. The partition function then becomes, 

z = £{v i( t) } / D{m(t)} J D{m} I{Dfa(t)} 



exp I fjj dt J2i iriiifa + 2vrt; 0i ) -^EijOi - fi jB )C , ii 1 (n J - - n x ) (44) 



- Ei 2§J + ipi ■ (Afa - A + 2^) 



where the path integral should be performed stepwize ||21|| . For the brevity of this expression 
we omitted the explicit time dependence of n», fa,Pi,Vi in the stepwize integrated action. 
This form allows us to understand the physical significance of p. The only ^-dependent 
term in the action is —ipi ■ A^. The derivative of the Lagrangian with respect to A is the 
current. Thus, pi is the Josephson current flowing through the site i. 

The path integral over the phase variables <pi(t) can now be performed. The phase 0, at 
the site i is coupled to the charge rii (via the term in^i) and to the vector p at the site i 
and its nearest neighbors. The integration over fa yields conservation of charge constraint 
on the integration over rii,Pi, in the form 

At-m + A ■ p t = (45) 

where the definition A t rii = — [rii(t + To) — rii(t)} makes the difference operator A a 3-vector. 

The constraint Q4"5p is nothing but a discretized form of a two dimensional conservation 
of charge equation. Like the latter, it can be solved by defining a 3-vector field K, that 
satisfies, 

e a ^A a K,^ = 2vrp JiM (46) 

where pi = (ni,Pi iX ,Pi >y ). The three components of /Q are real, like those of p^. The 
definition ( ^6] ) of /Q is not unique and it becomes unique only when a gauge is fixed. The 
partition function is, of course, independent of that gauge. The constrained path integral 
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over pi is replaced now by path integrals over /Q, constrained by the gauge condition. Here 
we choose to work in the Coulomb gauge, in which A ■ K, = 0. In that gauge the partition 
function becomes, 



Z = E Vi (t) / D{fC,} exp fg dt [-^Ey(A x jfyCjiK x K 



- E» s^jCA*^ + A/CgJ + < E,[A x (JC + /C«% • (v, - ^A,; 

(47) 

where K ext is defined by A x K, ext = 2nn x . We are now one step away from having an effective 
action for the vortices. The remaining step is an integration by parts of the last term in the 
action in (|4~7D . After performing that integration, the following action is obtained: 



S™ = $ dt J2i { i(j>r - + ilr ■ {Ki + K ext ) 

(48) 

+iE,(A x /Q)C^(A x Kj) + ^-((A,/Q) 2 + (A/C ,) 2 )} 

where the vortex 3-vector current J vor is defined as J vor = A x v, the average density of 
vortices is given by n v = £^ and the value of h has been restored. Equation (f4~8"D is the 
starting point of the discussion in section (3). 
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